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STATISTICS OF CENTROIDS OF VELOCITY 

A. Esquivel^ and A. Lazarian'^ 
RESUMEN 

Hacemos un resumen del uso de estadi'sticas de mapas de centroides de velocidad para obtener informacion de 
la turbulencia del medio interestelar a partir de observaciones. Los centroides de velocidad ban sido usados por 
mucho tiempo para obtener informacion de las propiedades de escalamiento de la velocidad turbulenta en el 
medio interestelar. Mostramos que, mientras que los mapas de centroides de velocidad son utiles para el estudio 
de turbulencia subsonica, no lo son para turbulencia supersonica, porque son influenciados por fluctuaciones 
de densidad. Mostramos tambien que para turbulencia sub-Alfvenica (tanto subsonica como supersonica) las 
estadi'sticas de dos puntos (como las funciones de correlacion o los espectro de potencias) muestran anisotropia. 
Esta anisotropia puede ser aprovechada para determinar la direccion del campo magnetico proyectado en el 
piano del cielo. 

ABSTRACT 

We review the use of velocity centroids statistics to recover information of interstellar turbulence from ob- 
servations. Velocity centroids have been used for a long time now to retrieve information about the scaling 
properties of the turbulent velocity field in the interstellar medium. We show that, while they are useful to 
study subsonic turbulence, they do not trace the statistics of velocity in supersonic turbulence, because they 
are highly influenced by fluctuations of density. We show also that for sub-Alfvenic turbulence (both supersonic 
and subsonic) two-point statistics (e.g. correlation functions or power-spectra) are anisotropic. This anisotropy 
can be used to determine the direction of the mean magnetic field projected in the plane of the sky. 
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1. INTRODUCTION 

It is well known that the interstellar medium 
(ISM) is turbulent. Such turbulence is magnetized 
and expands over several scales, ranging from au to 



kpc (Larson 1992 Armstrong, Rickett & Spangler 



1995; Deshpande, Dwarakanath & Goss 2000; Stan- 
[imirovic fc Lazarian 2001 Lazio et al. 2004). Under- 
standing of this magnetic turbulence is of great im- 
portance for key astrophysical processes, from star 
formation to diffusion of heat and cosmic rays, we 
refer the reader to recent reviews on the subject 



Elmegreen fc Scalo (20041, McKee fc Ostriker (2007). 



Observations of line-widths (Larson 1981 1992 



Scalo 


1984 




19871 and of the centroids of spectral 


lines 


1 von 


Hoerner| 1951 Miinch 


1958 Kleiner fc 


Dickman 1985 Dickman fc Kleiner 1985 Miesch fc 


Bally 


1994 


O'Dell fc Castaheda 


19871 have been 



used for well over half a century to study turbulence 
in the ISM. An important measure that one can hope 
to retrieve from spectral line data is the power-law 
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index of the underlying velocity field. However, the 
shape of spectral lines does not depend solely on the 
velocity field, but on the density of emitting mate- 
rial as well. The separation of the two contributions 
has proven to be a formidable (and in some aspects 
still an open) problem (see reviews by |Lazarian 2006 
^0081). 



In parallel with the increasing number and qual- 
ity of observations, there have been substantial theo- 
retical and numerical advances in our understanding 
of compressible magnetohydrodynamic (MHD) tur- 
bulence (see Cho fc Lazarian 2003 and also reviews 
Cho, Lazarian fc Vishniac (2003) and Cho fc Lazar- 
(2005)). Therefore comparison between theory, 



lan 



numerics and observations has become essential. 

Since turbulence is essentially an stochastic pro- 
cess, statistical methods are necessary for its study. 
Studies of correlations can be obtained in real space 
using correlation or structure functions, but also in 
Fourier space using spectra. Wavelets, which some- 
times are preferable for handling of the real observa- 
tional data, combine properties of both correlation 
functions and spectra and can be related to both of 
them. A well-known example of of wavelets, the A- 
variance, has been successfully used to retrieve the 
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same information yielded by two point statistics (i.e. 
power-spectra or structure functions) (Stutzki et al. 
T998: Mac Low fc Ossenkopf|2000|[Bensch et al.|2001| 



Ossenkopf & Mac Low 20021. However, neither of 



Coordinate Spectrum" (VCS; Lazarian & Pogosyan 
2006 2008). This paper explains when velocity cen- 



the measures above is capable of separating velocity 
and density contributions to the spectral lines. The 
latter separation should be done on the basis of the- 
oretical understanding of the statistical properties of 
the Doppler-shifted spectral linesj^ 

Recent years have been marked by a sharp in- 
crease of interest in statistical techniques of analy- 
sis of observations of astrophysical turbulence. We 
can mention in this respects two web sites by Alyssa 
Goodman, namely "Taste Tests"|2 where different 
comparisons of the numerical simulations and obser- 
vations are presented, and "Astronomical Medicine" 
sitc[^where application of sophisticated medical anal- 
ysis software are applied to astronomical data. In 
addition, we may mention new ways of analyzing 
column densities, such as "Genus" (Lazarian 1999, 
Lazarian, Pogosyan & Esquivel 2002, Kim & Park 
2007, Chepurnov et al. 2008), "Bispectrum" and 
"Bicoherence" (Lazarian 1999, Lazarian, Kowal & 
Beresnyak 2008, Burkhart et al. 2009). A discussion 
of those, is, however, beyond the scope of our present 
publication, which deals with a way of extracting the 
statistics of velocity from observed Doppler-shifted 
spectral lines. 

Among new theoretically-motivated way of re- 
covering velocity statists we can mention "Veloc- 



ity Channel Analysis" (VGA 


Lazarian & Pogosyan 


2000 


|Lazarian et al.||2001 


|2002| Esquivel et al. 


2003 


Lazarian &; Pogosyan 


2004 p "Modified Ve- 


locity Gentroids" (MVGs; Lazarian & Esquivel 


2003 


Esquivel & Lazarian||2005| Ossenkopf et al. 


2006 


Esquivel et al.||2007l the first three hereafter LE03 



and EL05, OELS06, respectively), and the "Velocity 

^An example of empirical approach to the problem is the 
use of the Principal Component Analysis (PCA) discussed, for 
instance, in Heyer & Brayn (2004). We feel that this powerful 
approach does not show all its strength being isolated from 
theory. For instance, it is known that for shallow density 
spectrum the statistics of Position-Position- Velocity (PPV) 
data cubes inseparably depends on both spectra of velocity 
and density (Lazarian &: Pogosyan 2000). This effect has not 
been demonstrated so far within the PCA approach. 

* www. cfa.harvard.edu/~agoodman/newweb/tastetests. html 

^am. iic.harvard.edu 

^The "Spectral Correlation Function" (SCF) l|Rosolowsky| 
|et al|1999[|Padoan, Rosolowsky &: Goodman|2001[ ) is another 
new measure, which, however, under close examination differ 
from the measure in VGA only by its normalization. The ad- 
vantage of the normalization adopted in the VCA is that the 
observations can be described in terms of underlying corre- 
lations of velocity and density, which is not the case for the 
normalization adopted for the SCF. 



troids, including MVGs are capable of measuring the 
statistical properties of the underlying velocity tur- 
bulence. 

The layout of this paper is as follows. In S|2]we 
present the basic statistical toolds used. We will re- 
view some of our work about the retrieval of velocity 
statistics from velocity centroids in ^ A special em- 
phasis on the anisotropies in the statistics that result 
from the presence of a magnetic field, and a discus- 
sion of how can they be used to reveal the direction 
of he mean magnetic field is presented in Sjlj Finally 
we provide with a brief summary in S|5] 

2. TWO-POINT STATISTICS 

Two-point statistics, such as correla- 
tion/structure functions and power spectra are 
the simplest, and most widely used method to char- 
acterize turbulence. The (second-order) correlation 
function of a quantity /(x) is defined as: 



5F(r) = ([/(x)-/(x + r)r 



(1) 



where r is the "lag", and (...) denotes an ensemble 
average over all the space (x). The correlation func- 
tion (or auto-correlation function) 



CF(r)-(/(x)./(x + r)), 



(2) 



can be easily related structure function as 
SF{r) = 2[CF{0)-CF{r)], where CF{0) is the 
variance of /(x). The power spectrum can be defined 
as the Fourier transform of the correlation function, 
P(k) = [CF{r)], where k is the wave-number. 

The usefulness of these type of functions lies in 
the fact that they have a power-law behavior in the 
so-called inertial range. Energy is injected at large 
scales, and cascades down without losses down to 
the scales at which (viscous) dissipation occurs. The 
inertial range is precisely between these two scales. 
For instance, the Kolmogorov model of hydrodynam- 
ical (and incompressible) turbulence predicts that 
the difference in velocities at different points in the 
fluid increases on average as the cubic root of the 
separation {\Sy \ cx l-^^^). This famous scaling results 
in a structure function SF{r) cx r^^^, and a (three 
dimensional) power spectrum P(fc) oc k~^^^^. Notice 
that the Kolmogorov model assumes isotropic turbu- 
lence (i.e. we have replaced r by r = |r|,and k by 
k — |k|). But, turbulence is not isptopic in general. 
In particular, by introducing a preferential direction 
of motion, the magnetic field that threads the ISM 
makes the turbulent cascade anisotropic. There are. 
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however, some statistical measures that are not very 
sensitive to such anisotropy (for instance VGA, see 
Esquivel et al.|2003 ). And in fact, it is customary to 
average correlation and structure functions in shells 
(or annuli) of equal separation (or average the power 
spectrum in wave number) effectively reducing the 
statistics to one dimension. In the following section 
(Sj3| we will do such averaging procedure, but we 
will come back to discuss how can these anisotropics 
be exploited to learn something about the magnetic 
field in g] 

3. TRACING THE STATISTICS OF VELOCITY 
WITH CENTROIDS 

Several studies, with varying degrees of success, 
have been made to obtain the spectral index of tur- 
bulence (power-law index of the velocity spectrum, 
correlation or structure function). However, many 
of them are restricted to ionized media (interestel- 
lar scintillations for instance), and more importantly 
they are only sensitive to fluctuations in density. 
While, density fluctuations are a natural result of 
a turbulent cascade, one has to make the leap from 
the observed fluctuations of density to a dynamical 
quantity predicted by theoretical models, such as ve- 
locity or magntetic field. It is therefore very desirable 
to obtain directly the spectral index of velocity from 
observations. 

Evidently, the Doppler-shifted spectral lines con- 
tain information about the velocity. However, they 
are also affected by density fluctuations, and one has 
to be careful to interpretate the statistics drawn from 
observations, in particular results from 2D maps of 
velocity centroids. Velocity centroids are ussually 
defined as ( von Hoerner||195T| Miinch|[l958 1: 



C(X) 



/ /line (X,?;^) dvz 



(3) 



where lu^e is the line intensity at a position X = 
{x, y) in the plane of the sky, at line of sight (LOS) 
velocity v^- The integration limits arc defined by 
the extent of velocities covered by the object. // 
the medium is optically thin, and the emissivity is 
proportional to the density (i.e. HI for instance), one 
can replace the velocity integrals by integrals over the 
actual LOS (chosen here to coincide with the z— axis, see 
LE03, EL05): 



C(X) 



/ (x) p (x) dz 
J P (x) dz 



(4) 



where x — {x,y,z). One could construct the struc- 
ture function of these centroids, but the denominator 
of eq. Q makes the algebra a bit messy, and does not 



provide a significant difference over "unnormalized cen- 
troids" (see LE03), 

S(X) = J Vz IlincOi.,V:,) dV:, ^ a J Vz{x.) p{ji) dz . (5) 

Replacing Xi = x and X2 = x + r, and similarly 
Xi — X and X2 = X + R, the structure function of 
centroids can be written as: 

{[S{X^) - S{X2)f) =0" jj dz,dz2 [d{v) - D{v)\^^^^^ 

(6) 

where 

D{V) = {[p (xi) (Xl) - p (X2) V. (X2)]') . (7) 

The notation jxi=X2 indicates that the integral is to be 
computed for a zero distance between Xi and X2 for 
the second term (i.e. varying only in z). Writing the 
density and velocity fields as a mean plus a fluctuating 
part (p = po + p, Vz — vo + Vz), and approximating 
the fourth order moments as a combination of second 
order moments (see LE03, EL05), the structure function 
of unnormalized centroids becomes: 

([^(Xi) - S{X2)f) ^ J1(R) + /2(R) + J3(R) + /4(R), 

(8) 

where 



J1(R) = a' {vl) Jfdz,dz2 [dp(r) - d,(r)\^^^^J^ ,(9a) 
J2(R) = <p2) JJdz^dz2 [d„,(r) - d,Ar)\^,=^,]{9h) 



73 (R) = 



/4(R) = 



1 „,2 
' 2 



// dzidz2 [dp{r)dv^ ( 



'^p('^)lxi=X2 ('■)lxi=X2 I ' 



(9c 



7/d2id^2 [c(r)- c(r)|^^^xj . (9d) 



We have made use of the 3D structure functions of the 
density and of the LOS velocity: 

dp(r) = ([p(xi)-p(x2)]'), (10a) 
dvA^) = ([S.(xi)-«. (X2)]'), (10b) 

and the remaining density-velocity cross-correlations 
have been grouped into 

c(r) = 2 {u3(xi)p(x2))^ - 4po {p(xi)?}2(xi)i;z(x2)) . 

(11) 

In the decomposition proposed in eq. (|8| the term II 
can be identified with the structure function of column 
density (weighted by («^), thus measurable from obser- 
vational data), the 12 term contains the analogous in 
terms of velocity, which could be used directly (assum- 
ing isotropy) to obtain the velocity spectral index. The 
remaining terms (13 and 14) are cross-terms, and cross- 
correlations, respectively, and they "contaminate" our 
statistics. Neglecting them is that we arrived to our def- 
inition of MVCs (LE03), we proposed to simply subtract 
the structure function of column density from that of the 



4 



ESQUIVEL & LAZARIAN 



centroids. Later on (EL05, OELS06,[E"squivel et al.|2007 1 



we tested the retrieval of the velocity spectral index with 
various data cubes, magnetohydrodynamic (MHD) sim- 
ulations and ensembles of artificially produced fractional 
Brownian motion cubes (fBms). Our results show that 
centroids are only useful to trace the scaling of turbu- 
lent velocity only when this is subsonic, and that the 
technique could be pushed to mildly supersonic turbu- 
lence (sonic Mach number < 2.5) with MVCs. More re- 
cently (Esquivel et al. 20071, we have found that the non- 



Gaussianity of both density and velocity fields (which is 
characteristic of highly supersonic turbulence) makes the 
approximation in eq. ([8| inadequate, the reason being 
that the fourth order correlations could not be well ap- 
proximated by second order ones in that case. As a guide- 
line for observers, in LE03 we proposed a necessary con- 
dition for centroids to trace the statistics of velocity (ne- 
glecting 13 and 14). When ([^(Xi) - ^(Xs)]^) /II < 1, 
then velocity centroids trace the velocity scaling (this is 
the case of subsonic turbulence), if the condition is only 
partially fulfilled (e.g. the ratio < 2) then MVCs would 
work, while the other centroids would not (this was the 
case of weakly supersonic turbulence simulations). The 
condition is strictly speaking a function of the lag R, 
thus one could have regions where it is true and regions 
in where it is not, depending on the slope of the structure 
functions involved. In fact one should compare them at 
the region in which we measure the spectral index. But, 
as a first approach, one can approximate the condition 
with the maximal value of the structure functions, lead- 
ing to {S^) / {P) <^ 1, which can be easily obtained from 
observations. A more restrictive criterion (but that re- 
quires additional observational information) is the ratio 
of the density dispersion to the mean density {jjp/po, 
OELS06), if this is less than unity centroids work, if not, 
then centroids are useless to get the velocity spectral in- 
dex. 

3.1. Application to SMC data 

To illustrate what we describe above, we have taken 
Ifl 21 cm. data that combines observations with the 
Parkes telescope and the Australia Telescope Compact 



Array (ATCA) interferometer, obtained by Stanimirovic 
et al. (19991. We obtained from them maps of column 



density, velocity centroids (normalized and unnormal- 
ized), and MVCs, they are shown in Figure [l] From 
both maps of centroids it is evident the global large scale 
motion of the SMC, a rotation respect to an axis that 
forms an angle of > 45° with respect to the x-axis in 
the maps of the Figure. On top with this regular motion 
there is an important turbulent velocity field. 

We have computed the structure function of these 
quantities, and we show them in the last panel of Fig- 
ure [l] For this object {S'^)/{P) ~ 0.36, thus we should 
expect centroids not to trace the scaling of velocity 
but rather should be highly affected by density fluc- 
tuations. This can be verified from the figure: II = 
(u^) ((Ji — I-z)^) (i.e. (n^) times the structure function 



of column density) is larger than any other quantity plot- 
ted, over the entire range of scales. The conclusion of 
this excercise is that to obtain the spectral index of tur- 
bulent velocity in the SMC one should resort to other 



techniques, such as VCA (see Stanimirovic & Lazarian 
200T|). 



4. THE EFFECT OF MAGNETIC FIELD: 
ANISOTROPY IN STRUCTURE AND 
CORRELATION FUNCTIONS 

The presence of a magnetic field introduces a prefer- 
ential direction of motion that makes the turbulent cas- 
cade anisotropic. In a magnetized plasma, the eddies 
become elongated in the direction of their local magnetic 
field (see for instance |Cho, Lazarian, fc Vishniac||2002| . 
Fortunately, in spite of being innadequate to recover the 
spectral index of velocity, velocity centroids refiect the 
anisotropic cascade in sub-Alfvenic turbulence (regard- 
less of the sonic Mach number). In isotropic turbulence, 
two-point statistics depend only in the magnitude of the 
lag (or wave- number) , therefore isocontours of 2D two- 
point statistics are circular. In magnetized turbulence, 
however, such isocontours become ellipses that are al- 
ligned with the mean magnetic field. The magnitude of 
the magnetic field determines how ellongated they are. 
Velocity statistics have been proposed before to study 
this technique (e.g. [Esquivel et al.|[2003| [Vestuto, Os-| 
triker fc Stone|2003| EL05), but to our knowledge, they 
have not been exploited observationally. 

In Figure [5] we show (taken from EL05) contours of 
equal correlation in maps of centroids from MHD sim- 
ulations with different sonic and Alfven Mach numbers, 
as well as different plasma (i (the ratio of gas to mag- 
netic pressures). The parameters of the simulations are 
included in the label of each panel, and in more detail in 
EL05. 

It is evident from the Figure the clear anisotropy 
present in all the sub-Alfvenic cases, regardless of the 
large differences in sonic Mach numbers and/or plasma 
p. The only case in which the anosotropy is not evident is 
for super- Alfvenic turbulence (panel [c], with Ma ~ 8). 
Whether turbulence in the ISM is typically sub-Alfvenic 
or super- Alfvenic is an open debate, for instance [Padoan] 
et al.| (j2004 ) advocate for a model of supersonic turbu- 



lence in molecular clouds. 

A different way to visualize the anisotropy in two- 
point statistics is to separate the correlations into per- 
pendicular and parallel to the mean (i.e. global) mag- 
netic field. This would mean to plot the value of the 
correlations only along the symmetry axes of the ellipses 
of equal correlation. We present this in Figure |3] (from 
EL05). 

In this case the anisotropy reveals itself as two dis- 
tinct correlation lenghts for the parallel and perpendic- 
ular direction. The difference in scale-lengths should in 
principle refiect a dependence on (B/Bq)^, where B is 
the fluctuating magnetic field and Bq the mean mag- 
netic field. One should have in mind, however, that this 
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Fig. 1. Application to Small Magellanic Cloud observations (data from Stanimirovic et al. 1999| ) Top left: column 
density map. Top right and bottom left: unnormalized velocity centroids and normalized velocity centroids, respectively. 
Bottom right: structure functions of the maps of column density and centroids (as indicated in the legend). 



method is only sensitive to the direction of the magnetic 
field in the plane of the sky. Turbulence with a strong 
mean magnetic field oriented along the line of sight would 
be indistinguishable from super-Alfvenic turbulence. It 
is therefore difficult to determine (B/Bq) from this type 
of statistics in supersonic turbulence (where density fluc- 
tations are important). Nonetheless, the results in terms 
of the direction of the Bo are robust and could be used 
when other methods are not available. 



5. DISCUSSION AND SUMMARY 

We have made a review of our previous work on veloc- 
ity centroids, and their application to retrieve the scaling 
properties of the turbulent velocity field (i.e. spectral in- 
dex). Additional details can be found in Esquivel et al. 
(20031); LE03; EL05; OELS06, and Esquivel ct al. (2007j). 



This work on centroids assumes a optically thin medium 
with emissivity proportional to density (e.g. HI), while 
self absorption has been addressed within VGA (Lazar- 



ian & Pogosyan 2004 1 and VCS (Lazarian & Pogosyan 



2006) techniques, which opens a way of formulating the 
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a) M3 ~ 0.5, Ma ~ 0.7, |S = 4 b) M3 ~ 2.5, ~ 0.7, ^ = 0.2 
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Fig. 2. Taken from EL05: Anisotropy in the correlation functions: contours of equal correlation for MHD simulations 
{solid lines). For reference we show isotropic contours as dotted lines. The sonic and Alfven Mach numbers {A4s, Ma 
respectively), and the plasma arc indicated in the title of the plots. The anisotropy reveals the direction of the 
magnetic field for all the sub-Alfvenic cases, regardless of the sonic Mach number. 



theory of centroids for partially absorbing media. This 
has not been done yet. 

We also assumed that the emission lines are being 
used. If the turbulence volume is between the observer 
and an extended emission source, centroids can be also 
used with the absorption lines. The present theory as- 
sumes that the absorption is in linear regime however. 
Needless to say, that extend of the emission source of 
multiple emission sources will determine the spatial cov- 
erage of scales that are testable with the centroids. In 



comparison, the VCS technique can deal with saturated 
absorption lines and does not depend on the spatial cov- 
erage of the data. 

The advantage of centroids compared to the VGA 
and the VCS techniques is, first of all, the ability of 
centroids to study magnetic field direction, and, second, 
ability to study subsonic turbulence more reliably. While 
the VGA and VGS are also capable of studying subsonic 
turbulence, the procedures for extracting of information 
within these techniques are much more complicated. 
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Fig. 3. Taken from EL05: Correlation functions taken in directions parallel and perpendicular to the mean magnetic 
field. The anisotropy shows in the different scale-lengths for the distinct directions. It is noticeable the little difference 
in panels (a), (b), and (d). Which correspond to the same ratio B/Bq, but very different sonic Mach numbers. Panel 
(c) corresponds to super-Alfvenic {B > Bo) turbulence, and the anisotropy is not evident in the centroid maps. 



The most frequent mistake, we feel, is the use of cen- 
troids while dealing with supersonic data (see Miville- 
Deschenes et al. 2003). It is important to under- 
stand that when wo deal with the multi-phase interstellar 
medium, e.g. HI, the criterion of being subsonic is the 
most difficult to be satisfied for the cold medium. There- 
fore, our research shows that while the application of ccn- 
troids to fill regions is justified, their applications to find 
the statistics of velocities in molecular clouds and cold 
HI cannot deliver reliable spectra. This, as we shown, 
should not discourage the use of the centroids for stud- 
ies the direction of magnetic field. The reliability of the 
latter technique can be tested by comparing the polar- 



ization arising from aligned dust (see Lazarian & Hoang 
2007 for a discussion when we can rely on grain alignment 
to trace magnetic fields) and the results of the statistical 

anisotropy analysis with centroids. 

The main results of the paper above can be summa- 
rized as the following: 

• Centroids maps can be used to trace the spectral 
index of the underyling turbulent velocity in sub- 
sonic turbulence. In mildly supersonic turbulence 
(sonic Mach number < 2 can be studied with MVCs 
(LE03, EL05). 

• Two criteria can be used to determine if centroids 
are useful to obtain the velocity spectral index. If 
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([5(Xi) - ^(Xa)]^) //I < 1 (see eqs 8, 9) or if 
OpI pa < 1. The first is a necessary condition, the 
latter is more robust measure, but might require ad- 
ditional information than what is ussually available 
from spectroscopic observations. 

We presented an example of the application of veloc- 
ity centroids to try retrieving the velocity spectral 



index from real data (SMC observations from Sta- 
[nimirovic et al.|1999[ ). The criteria proposed above 
was not fulfilled, thus the spectral index should be 
studied using an alternate method (such as VGA, 
see |Stanimirovic fc Lazarian|200l| . We showed how 
the dominant contribution to the statistics of cen- 
troids are density fluctuations. 

• Two point statistics are anisotropic for sub-Alfvenic 
turbulence, both for subsonic and supersonic cases. 
This anisotropy is evident if we plot isocontours of 
velocity centroids, which become elongated and re- 
veal the direction of the mean magnetic field pro- 
jected onto the plane of the sky. This anisotropy 
can be used where other magnetic field measures 
are not available. 
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